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Various  properties  of  vitrifying  liquids  are  correlated  with  the  dispersity  of  the  dynamics,  the  latter 
reflected  in  the  magnitude  of  the  nonexponentiality  parameter,  (3K ,  describing  the  distribution  of 
relaxation  times.  These  properties  include  the  mean  relaxation  time,  ra,  the  fragility,  and  the 
dynamic  crossover.  The  correlations  with  f3K  are  observed  in  both  experimental  data  and  the  results 
from  molecular  dynamics  simulations  on  Lennard-Jones  (LJ)  type  systems.  Another,  rather  obvious 
property  to  correlate  with  j3K  is  the  dynamic  heterogeneity,  which  can  be  quantified  from  the  number 
of  molecules,  Nc,  dynamically  correlated  over  a  time  span  ra.  For  a  given  LJ  system,  Nc  can  be 
rigorously  calculated  and  we  find  that  it  does  indeed  correlate  with  /3K  over  a  range  of 
thermodynamic  conditions.  However,  the  analysis  of  experimental  data  for  a  broad  range  of  real 
materials,  wherein  an  approximation  is  required  to  obtain  Nc,  reveals  the  absence  of  any  relationship 
between  Nc  and  (iK  among  different  materials.  ©  2010  American  Institute  of  Physics. 
[doi:  10. 1063/1. 348 1355] 


I.  INTRODUCTION 

Viscous  liquids  are  sufficiently  dense  for  molecules  to 
exert  a  reciprocal  influence,  and  these  correlated,  many-body 
interactions  distinguish  the  supercooled  regime  from  the 
much  simpler  dynamics  at  high  temperature.  The  degree  of 
intermolecular  cooperativity  varies  spatially,  with  the  dy¬ 
namic  correlations  persisting  over  time  scales  on  the  order  of 
the  relaxation  time.  Since  dynamic  heterogeneity  is  intrinsic 
to  the  supercooled  liquid  state,  it  is  an  obvious  quantity  to 
characterize  dynamic  properties. 1-4  One  consequence  of  dy¬ 
namic  heterogeneity  is  deviation  of  the  linear  susceptibility 
from  the  Debye  or  exponential  behavior,  with  the  relaxation 
function  becoming  more  dispersive  to  reflect  the  heterogene¬ 
ity  of  the  local  dynamic  environments.  A  variety  of  experi¬ 
mental  methods  have  been  employed  to  demonstrate  that 
nonexponential  relaxation  results  from  a  superposition  of  dy¬ 
namically  distinguishable  contributions.5  This  nonexponenti- 
ality,  commonly  assessed  from  the  magnitude  of  the  Kohl- 
rausch  stretch  exponent,  (3K ,  is  known  to  correlate  with 
relaxation  properties  such  as  fragility6  and  the  dynamic 
crossover  (change  in  dynamics  above  Tf.7  Additionally,  [3K 
is  related  to  phenomena  such  as  rotational-translational 
decoupling  and  the  confinement-induced  enhancement  of 
molecular  mobilities.  1--14  Thus,  quantifying  the  dynamic  het¬ 
erogeneity  and  nonexponentiality  of  liquids  can  provide  in¬ 
sight  into  the  fundamental  variables  governing  the  dynamics 
near  Tg.  Indeed,  models  of  the  glass  transition  generally  posit 
a  growing  dynamic  correlation  length  as  causing  the  marked 
increase  of  ra  in  vitrifying  liquids15-19  (although  recent  simu¬ 
lation  results  suggest  an  alternative  viewpoint-0). 
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Since  the  conventional  description  of  the  dynamics  in¬ 
volves  the  correlation  of  a  variable  at  two  times  (for  ex¬ 
ample,  ( p(k,0)p(k,t )),  where  p  is  the  density  and  k  is  the 
wave-vector),  it  follows  that  a  four-point,  time-dependent 
correlation  function  is  required  to  describe  dynamic 
heterogeneities.”  ’  Thus,  the  susceptibility  defined  in  terms 
of  spatial  and  temporal  correlations, 

X4(0  =  J  (p(r i>0)p(r ,  +  r2,0)p(rl,t)p(rl  +  r2,t))rldr2, 

(1) 


has  a  maximum  at  t  ~  ra  that  is  proportional  to  Nc,  the  num- 

22—24 

ber  of  correlating  molecules-” 

Vf  =  max{x4«}  =  *r.  (2) 


Efforts  to  study  dynamic  heterogeneities  have  been  stymied 
by  the  difficulty  of  determining  X4 W  since  conventional  re¬ 
laxation  spectroscopies  yield  the  linear  susceptibility.  How- 
ever,  molecular  dynamics  simulations”  ’  and  multidi¬ 
mensional  NMR”9  12  have  been  employed  to  obtain  ^4(f) 
respectively  for  model  particle  systems  and  a  few  real  mate¬ 
rials  (data  for  the  latter  limited  to  low  temperatures). 

A  development  in  quantifying  dynamic  heterogeneities 
in  real  materials  was  the  derivation  by  Berthier  et  al.  of 
xft)  in  terms  of  the  temperature  derivative  of  two-point  lin¬ 
ear  susceptibilities,  <J>(f), 


*4  0) 


dT  j  ’ 


(3) 


where  A cP  is  the  isobaric  heat  capacity  change  at  T „  and  kB 
is  the  Boltzmann  constant.  This  offers  the  possibility  of  ob- 


0021-9606/201  0/1 33(1 2)/124507/4/$30.00  133,  124507-1 


©  2010  American  Institute  of  Physics 


Downloaded  29  Oct  2010  to  132.250.22.10.  Redistribution  subject  to  AIP  license  or  copyright;  see  http://jcp.aip.org/about/rights_and_permissions 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

2010 

2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2010  to  00-00-2010 

4.  TITLE  AND  SUBTITLE 

Correlation  Of  Nonexponentiality  With  Dynamic  Heterogeneity  From 
Four-Point  Dynamic  Susceptibility  X4(t)  And  Its  Approximation  XT(t) 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory, Washington, DC, 20375 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

The  Journal  Of  Chemical  Physics  133,  2010,  Government  or  Federal  Purpose  Rights  License. 


14.  ABSTRACT 

Various  properties  of  vitrifying  liquids  are  correlated  with  the  dispersity  of  the  dynamics,  the  latter 
reflected  in  the  magnitude  of  the  nonexponentiality  parameter,  K,  describing  the  distribution  of  relaxation 
times.  These  properties  include  the  mean  relaxation  time,  ,  the  fragility,  and  the  dynamic  crossover.  The 
correlations  with  K  are  observed  in  both  experimental  data  and  the  results  from  molecular  dynamics 
simulations  on  Lennard-Jones  LJ  type  systems.  Another,  rather  obvious  property  to  correlate  with  K  is 
the  dynamic  heterogeneity,  which  can  be  quantified  from  the  number  of  molecules,  Nc,  dynamically 
correlated  over  a  time  span.  For  a  given  LJ  system,  Nc  can  be  rigorously  calculated  and  we  find  that  it  does 
indeed  correlate  with  K  over  a  range  of  thermodynamic  conditions.  However,  the  analysis  of  experimental 
data  for  a  broad  range  of  real  materials,  wherein  an  approximation  is  required  to  obtain  Nc,  reveals  the 
absence  of  any  relationship  between  Nc  and  K  among  different  materials. 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

ABSTRACT 

18.  NUMBER 

OF  PAGES 

19a.  NAME  OF 

RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

5 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


124507-2 


Roland  et  al. 


J.  Chem.  Phys.  133,  124507  (2010) 


taining  the  number  of  molecules  with  correlated  dynamics  in 
terms  of  readily  accessible  experimental  quantities, 

NC(T,P)  ~  T^T^max  Xi(t,T,P)}2.  (4) 

A  cP(P)  i 

A  number  of  works  have  appeared  in  which  the  XtU)  aP~ 
proximation  was  used  to  evaluate  Nc  for  different 
materials'  ’  and  different  thermodynamic  conditions.  A 
comparison  of  simulation  results  for  xa  and  Xt  indicated  rea¬ 
sonable  correspondence  between  the  two  methods  for  longer 
27 
v 

Conceptually  the  spatial  extent  and  amplitude  of  the  dy¬ 
namic  correlations  are  expected  to  have  a  connection  to  the 
dispersion  in  the  relaxation  times  since  dynamic  heterogene¬ 
ity  is  defined  by  molecular  mobilities.  By  the  same  token 
properties  that  correlate  with  f3K  should  also  correlate  with 
Nc  and  vice  versa.  Up  until  now,  any  relation  between  f3K  and 
Nc  has  not  been  evaluated.  In  this  paper,  we  provide  the  test, 
using  simulations  of  the  Kob-Andersen  Lennard-Jones  (LJ) 
mixture  '6  as  a  function  of  temperature  T  and  density  p.  The 
nonexponentiality  parameter  [3K  is  obtained  from  the  self¬ 
intermediate  scattering  function  and  Nc  from  the  four-point 
dynamic  susceptibility,  XaU)-  Correlation  between  (3K  and  Nc 
is  assessed  by  comparing  the  changes  of  these  two  quantities 
on  varying  the  thermodynamic  conditions,  i.e.,  at  different 
combinations  of  T  and  p.  To  test  the  correlation  between  [3K 
and  Nc  obtained  from  Xi  via  Eq.  (4),  we  use  experimental 
data  of  45  glass-formers  including  polymers,  oxide  glass- 
formers,  selenium,  hydrogen  bonded  materials,  and  van  der 
Waals  glass-formers. 

II.  RESULTS 

A.  Correlation  between  stretching  parameter  fiK  and 
Nc  from  x?a* 

28 

Previously  it  was  shown”  from  molecular  dynamics 
simulation  of  the  Kob-Andersen  binary  mixture  with  par¬ 
ticles  interacting  according  to  a  12-6  LJ  potential  that  the  full 
f-dependence,  and  thus  the  maximum  in  /y4(f)(=A/r),  were 
invariant  over  a  relevant  range  of  T  and  p  for  state  points  for 
which  the  scaling  variable  pyIT  is  constant  with  7=5.07. 
Moreover,  the  reduced  relaxation  time  t*,  defined  as  the 
Kohlrausch  decay  time  for  the  self-intermediate  scattering 
function  Fs(k,t )  multiplied  by  p1/37’1/2,37  is  also  invariant  for 
state  points  for  which  pr/  T  is  constant.  Thus,  the  same  value 
of  the  material  constant  y=5.07  superposes  both  y4ax  and  t 
of  the  system  versus  pyl  T,  which  means  that  the  number  of 
dynamically  correlated  particles  is  a  unique  function  of  r*,  at 
least  in  the  regime  where  density  scaling  holds. 

Values  of  /3K  were  obtained  by  nonlinear  fitting  of  the 
stretched  exponential  function  <J>(r)  =  «I>0  exp(-(f/  t)@k)  to  the 
long-time  portion  of  Fs(k,t)  at  a  reduced  wave-vector  kp1/3 
=  7.44  over  the  same  range  of  T  and  p  as  in  Ref.  28.  The  fit 
was  cut  off  at  short  times  where  the  contribution  of  the  fast 
initial  decay  of  Fs(k,t)  became  significant;  that  is,  where  the 
shape  of  the  curve  begins  to  deviate  from  a  stretched 
exponential.  Representative  fits  are  shown  in  the  inset  of 
Fig.  1.  This  figure  shows  f3K  plotted  as  a  function  of  pyIT 
using  y=5.07,  the  same  exponent  superposing  versus 


FIG.  1.  Kohlrausch  exponent  as  a  function  of  the  product  variable  p5  01/T 
for  simulated  LJ  particles.  Each  symbol  represents  a  distinct  state  point.  Fits 
to  the  self-intermediate  scattering  function  at  longer  time  are  shown  for  two 
temperatures  in  the  inset.  Simulation  data  are  from  Ref.  28. 

pyIT  in  Ref.  28.  Within  the  experimental  error,  (3K  is  invari¬ 
ant  for  state  points  for  which  the  scaling  variable  pyIT  is 
constant.  In  Fig.  2  these  same  (3K  are  plotted  as  a  function  of 
t*,  and  the  results  indicate  that  [3K  is  a  unique  function  of  r* 
for  any  thermodynamic  condition  in  the  regime  where  den¬ 
sity  scaling  holds.  This  constancy  of  f3K  for  all  state  points  at 
any  fixed  r„  has  been  shown  from  experimental  data  on 
many  glass-formers.39'40 

To  compare  dynamic  correlation  to  nonexponentiality,  in 
Fig.  3  y4lax  is  plotted  versus  (3K.  The  correspondence  is  good 
and  in  accord  with  the  usual  interpretation  of  dynamic 
heterogeneity — an  increasing  correlation  (larger  Nc)  associ¬ 
ated  with  an  increasing  breadth  of  the  relaxation  function 
(smaller  (3K).  Also  shown  in  Fig.  3  are  the  values  of  Nc 
calculated  from  Xii1)  [Eq.  (4)].  The  agreement  is  good  at  low 
temperatures  and  higher  densities,  but  otherwise  the  approxi¬ 
mation  underestimates  the  number  of  dynamically  correlated 
particles.  This  agrees  with  previous  simulation  results  com¬ 
paring  XaU)  and  Xii1)  as  a  function  of  temperature  at  zero 
27 

pressure. 


FIG.  2.  Kohlrausch  exponent  as  a  function  of  the  relaxation  time  determined 
from  fits  of  Fs(k,t)  at  longer  time.  Each  symbol  represents  a  distinct  state 
point.  Simulation  data  are  from  Ref.  28. 
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FIG.  3.  Nc  from  Xa  [Eq.  (2) — open  symbols]  and  from  Xt  [Eq.  (4) — filled 
symbols]  vs  the  Kohlrausch  exponent  for  LJ  particles.  Each  symbol  repre¬ 
sents  a  distinct  state  point.  Increasing  dynamic  heterogeneity  is  associated 
with  a  broader  dispersion  of  relaxation  times.  The  Xt  approximation  departs 
from  the  rigorous  calculation  at  higher  temperatures  and  lower  densities. 
Simulation  data  are  from  Ref.  28. 


B.  Lack  of  correlation  between  [iK  and  Nc  from  *7 (f) 
approximation 

In  Fig.  4  plotted  versus  fiK  is  Nc  calculated  at  Tg  from 
Eq.  (4)  using  experimental  measurements  for  45  materials. 
The  Nc  data  were  taken  from  Ref.  34  and  (iK  are  as  reported 
in  the  original  papers  (references  in  Ref.  34).  Since  a  larger 
number  of  dynamically  correlated  molecules  should  coincide 
with  a  greater  dispersity  of  relaxation  times,  we  expect  these 
quantities  to  be  inversely  related.  However,  a  linear  correla¬ 
tion  of  the  data  in  Fig.  4  yields  a  Pearson  correlation  coeffi¬ 
cient  R  =  —  0.22,  indicating  an  absence  of  dependence.  If  we 
replace  Nc  by  the  corresponding  molar  volume  or  the  number 
of  chemical  groups  (e.g.,  “beads”  as  defined  by  Stevenson 
and  Wolynes41),  the  correlation  deteriorates  further,  with  the 
absolute  value  of  R  becoming  smaller. 

If  <l>(f)  has  the  Kohlrausch  form,  to  a  very  good  approxi¬ 
mation  Eq.  (4)  can  be  rewritten  so  that  (3K  appears 
explicitly,34 


FIG.  4.  Plot  of  Kohlrausch  exponent  vs  Nc  obtained  from  the  Xi(0  approxi¬ 
mation  for  45  glass-formers  at  r?  and  atmospheric  pressure  (data  are  from 
Ref.  34):  polymers  (diamonds),  oxide  glass-formers  and  selenium  (circles), 
hydrogen  bonded  materials  (triangles),  and  van  der  Waals  glass-formers 
(squares).  The  absence  of  correlation  is  indicated  by  the  small  value  of  the 
Pearson  linear  correlation  coefficient. 


FIG.  5.  Fragility  corresponding  to  data  in  Fig.  4.  Correlation  coefficient  is 
indicated. 


b  (  PkYI  d  In  ra\ 

A  cP  \  e  )  \  d  In  T  /  ’ 


where  e  is  Euler’s  number.  One  of  the  most  common  dy¬ 
namic  properties  of  supercooled  liquids  is  their  fragility,  m, 
proportional  to  the  apparent  activation  enthalpy  at  Tg  normal¬ 
ized  by  kBT\  that  is,  m  is  proportional  to  the  last  factor  in 
brackets  in  Eq.  (5).  It  is  well  accepted  that  m  correlates  with 
the  breadth  of  the  relaxation  time  distribution;  this  has 
even  been  quantified  as  ;h  =  250(±  30) -320/3^  (Ref.  6)  or 
in'* {$£  (Ref.  34).  This  implies  that  a  smaller  f3K  in  Eq.  (5)  is 
compensated  by  a  larger  fragility.  The  inference  is  that  the 
variation  of  Nc  among  different  materials  is  almost  entirely  a 
consequence  of  differences  in  A cP.  Thus,  both  the  form  of 
Eq.  (5)  and  the  data  in  Fig.  4  indicate  that  there  is  essentially 
no  relationship  between  the  dispersion  of  relaxation  times 
and  the  degree  of  correlation  of  the  dynamics  as  quantified 
using  Xiif )■ 

Since  /3K  and  m  are  correlated  (inversely),  it  follows 
from  the  results  in  Fig.  4  that  Nc  and  m  would  exhibit  no 
mutual  dependence.  Previously  Berthier  et  al.  examined 
data  for  15  materials  and  concluded:  “Dynamic  correlations 
revealed  by  Xiif)  increase  weakly  with  fragility.”  Extending 
the  data  set  threefold  (Fig.  5)  emphasizes  the  weakness  of 
this  putative  correlation.  Hong  et  al.4~  recently  reported  that 
the  dynamic  heterogeneity  length  scale  deduced  from  the 
frequency  of  the  Boson  peak  for  various  polymers  and  mo¬ 
lecular  liquids  is  not  correlated  with  fragility.  This  is  consis¬ 
tent  with  our  results  herein.  Moreover,  at  Tg  Nc  from  Eq.  (5) 
shows  an  inverse  proportionality  to  the  configurational  en¬ 
tropy  (assuming  the  latter  is  related  to  the  heat  capacity 
change  at  7’,,), 3 1  which  implies  no  correlation  with  fiK. 


III.  CONCLUSIONS 

We  find  that  max{^4(f)}  correlates  with  (iK  as  T  and  p  are 
varied  for  a  given  system,  and  although  the  Xt  approximation 
is  inaccurate  away  from  Tg,  using  the  latter  to  obtain  Nc  this 
correlation  with  (3K  is  maintained  for  a  given  material.  The 
problem  arises  when  comparing  different  systems  at  their 
respective  Tg.  Different  materials  exhibit  different  Nc  (from 
Xt)  for  a  given  fiK  because  of  differences  in  Ac,,.  The  form 
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of  Eq.  (4)  obviates  the  possibility  of  a  correlation  of  dynamic 
heterogeneity  with  either  /3K  or  fragility. 

Thus,  it  appears  that  short  of  abandoning  the  idea  that  a 
growing  xa  susceptibility  implies  more  dispersive  dynamics, 
it  may  be  necessary  to  re-examine  the  assumptions  underly¬ 
ing  Xt-  Its  derivation  assumes  that  the  contribution  to  the 
nonlinear  susceptibility  from  density  fluctuations  is  signifi- 
candy  smaller  than  the  effect  of  energy  fluctuations."  How¬ 
ever,  for  van  der  Waals  molecular  liquids,  density  and  ther¬ 
mal  energy  exert  about  equal  effects  on  the  linear  response, 
and  even  for  polymers  the  contribution  from  density  is  not 
negligible.43  Note  that  differences  between  the  exact  and  ap¬ 
proximation  values  for  Nc  are  larger  for  an  LJ  system  than 
for  simulated  silica  (at  least  at  higher  T),'1 2 * 4 5 6 7  and  indeed  the 
density  contribution  is  expected  to  be  more  important  for  the 
liquid  than  for  the  network  glass-former. 

Berthier  et  al.  pointed  out  that  the  susceptibility  can 
change  not  only  due  to  a  growing  length  scale  but  also  be¬ 
cause  of  the  growth  in  the  magnitude  of  the  local  fluctua¬ 
tions.  Whereas  f3K  depends  on  how  much  the  individual  (mo¬ 
lecular)  time  varies  within  this  volume,  greater  differences 
among  reorientational  rates  do  not  necessarily  imply  a  longer 
length  scale  for  decay  of  mutual  correlations.  In  other  words, 
how  different  two  quantities  are  is  not  equivalent  to  how  fast 
their  correlation  decays.  Nevertheless,  it  is  not  obvious  how 

the  amplitude  of  the  fluctuations  would  not  be  manifested  in 
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the  magnitude  of  [3K,  although  in  the  model  of  Wolynes, 
different  glass-formers  are  distinguished  by  the  range  of  ac¬ 
tivation  barriers  within  their  correlation  volume,  rather  than 
by  the  size  of  the  latter.  Clearly,  further  work  on  real  mate¬ 
rials  measured  by  techniques  such  as  multidimensional  NMR 
is  required  in  order  to  better  assess  the  relationship  between 
the  relaxation  time  dispersion  and  Nc,  and  the  utility  of  the 
Xt({)  approximation  in  quantifying  the  latter. 
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